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Abstract 

We study motion of a relativistic charged particle in a plane slow electromag- 
netic wave and background uniform magnetic field. The wave propagates 
normally to the background field. Under certain conditions, the resonance 
between the wave and the Larmor motion of the particle is possible. Capture 
into this resonance results in acceleration of the particle along the wave front 
(surfatron acceleration). We analyse the phenomenon of capture and show 
that a captured particle never leaves the resonance and its energy infinitely 
grows. Scattering on the resonance is also studied. We find that this scat- 
tering results in diffusive growth of the particle energy. Finally, we estimate 
energy losses due to radiation by an accelerated particle. 

1 Introduction 

Motion of a charged particle in an inhomogeneous electromagnetic field can 
be described by a Hamiltonian nonlinear system, which cannot be solved 
analytically in a general case [JJ. In such a system various resonant effects 
take a place [21 [3] , and some of them can by described using the adiabatic 
approach. In particular, a charged particle in the field of an electromagnetic 
(or electrostatic) wave and a background magnetic field can be trapped in the 
potential well of the wave and accelerated along the wave front. This phe- 
nomenon is called a surfatron acceleration. The mechanism of the surfatron 
acceleration of charged particles is often considered for description of various 
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plasma-physics phenomena [H [5]. Originally this mechanism was suggested 
for description of charged particles acceleration along the front of a shock 
wave jl] and this application is still actual [6]. On the other hand, there are 
various astrophysical applications of the surfatron acceleration mechanism 
to problems of generation of high energy particles [7J El [10] and conse- 
quent radiation [HI [HI 13] ■ Surfatron acceleration of relativistic particles 
was considered, for example, in JTUdS]. In all these papers authors consider 
a particle interaction with an electrostatic wave. 

Surfatron acceleration of a particle by an electromagnetic wave is less 
studied. The analytical theory was constructed only for nonrelativistic [16] 
or ultrarelativistic [HI [T7] particles. The effect of large particle velocity 
was estimated numerically in [HI [T7J. Also, numerical calculations were 
carried out for the case when the wave amplitude is small compared to the 
background magnetic field [HI [20] ; in this case, the particle is accelerated by 
multiple scatterings on the wave. In addition, several laboratory experiments 
with relativistic particles and large wave amplitudes were carried out for the 
investigation of surfatron acceleration of charged particles by electromagnetic 
waves (see [21] and references therein). Therefore, a complete analytic theory 
of relativistic charged particle captures by electromagnetic waves and the 
resulting acceleration is important. 

Particle capture and surfatron acceleration is possible if phase velocity of 
the wave is smaller than the absolute value of the particle velocity (and hence, 
smaller than the speed of light). In this case the projection of particle velocity 
onto the wave vector direction can become equal to value of the phase velocity 
of the wave and the resonance takes a place. There are several plasma modes 
that can support a wave with needed properties: the magnetosonic wave 
with frequency close to lower- hybrid [22J, the plasma wave with frequency 
close to higher-hybrid [7J or various drift modes of plasma instability [23J. 
In addition, a relatively small population of trapped particles can decrease 
the phase velocity of a wave [21] and establish the condition of resonance 
interaction. 

A secondary effect of surfatron acceleration is the radiation of acceler- 
ated particles due to the oscillatory component of their motion [HI |23]. A 
captured particle accelerates along the wave front, and at the same time it 
oscillates near the minimum of the wave potential well. Due to these oscil- 
lations the particle can radiate. Estimates of this radiation were carried out 
in the case of an electrostatic wave. The question of particle radiation in the 
system with electromagnetic wave is discussed in our paper. 
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We study the problem of interaction of a charged particle in a uniform 
magnetic field with an electromagnetic wave using the theory of resonant 
phenomena. The study of slow passages of a Hamiltonian system through a 
nonlinear resonance was started in [2] . In the present paper we use the theory 
of resonant processes in Hamiltonian systems with slow and fast motions in 
the form developed in E3 EE] (see [26J for references to preceding works) . 
The description of scattering on resonances and captures into resonances 
plays a central role in this theory. Resonant phenomena arise in a variety 
of problems of physics, including hydrodynamics, celestial mechanics, and 
plasma physics. For several examples of resonant phenomena, see, e.g., [29| 



2 Main equations 

We consider motion of a relativistic charged particle of mass m and charge 
e in a uniform magnetic field B = (0, 0, B ) and the field of a plane linearly 
polarized electromagnetic wave propagating in perpendicular direction to 
B . Thus, in the Cartesian coordinates (<?i, <?2, (fa) the resulting magnetic 
field components are B\ = B 2 = 0, B 3 = —B w sin(A;g 1 — ut) + B , where B w 
is the amplitude of the magnetic field of the wave, k is the magnitude of the 
wave vector directed along the gi-axis, and Cj is the wave frequency. The 
corresponding vector potential can be chosen as 

A = ^0, B oqi + ^j- cos(A;gi - cut), 0^ . (1) 

Let pi, i = 1, 2, 3 be components of the particle's momentum. Introduce 

e eB w 

V2 = P2 + --Bo9i H — cos(A;gi - tit). (2) 

c ck 

The Hamiltonian of the system is: 

H = ^ m 2 c 4 + c 2 pl + ^cV 2 — e£?o<?i — — ^ cos(/cgi — tui)^ + c 2 p\, (3) 

and pairs of canonically conjugate variables are (pi, qi), (P 2) ^2), (P3, 93)- The 
Hamiltonian does not contain variables g 2 and q%. Thus canonically conjugate 
momenta V2 and p 3 are constants of motion. We can put p 3 = (it always can 
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be done by redefining the particle's mass) ; one can also make V2 = choosing 
properly the origin in q x . Introduce Larmor frequency ul = eB /(mc) and 
dimensionless parameter e = eB w /(mkc 2 ). We assume that e is small: < 
e <^ 1. Use the following rescaling to make the system dimensionless: 

Pi u L qi 

Pi = , ft = , 

mc ec 

e mc A 

k CS CjE 
K = , U = . 

The Hamiltonian in the new variables is: 



H = + (eq x + e cos(A;g 1 - oot)) 2 . (4) 

Introduce new variable U = out. Let Pu be the variable, canonically conjugate 
to U. Thus we obtain a 2 d.o.f. Hamiltonian system. The Hamiltonian takes 
the form: 

U = uP u + sjl+p\ + (eqi + e cos(k qi - U)) 2 . (5) 

Now we introduce the wave phase ip = kq\ — U as a new variable. To this end, 
we make canonical transform with generating function W = pq± + I(kqi — U), 
where J is a new variable canonically conjugate to ip. For the new variables 
p, q, I we have q — qi, p — Pi — Ik, I = —Pu- Denote eq = q. Omitting the 
tilde, we find for the Hamiltonian in the new variables: 

H = -ujI+ y/l + (p + Ik) 2 + (q + e cos ^) 2 , (6) 

where the pairs of canonically conjugate variables are (I,<f) and (p, e~ 1 g). 
This Hamiltonian can be represented in the form H = H + eH\, where 

H = -ul + V 1 + (P + Ik ) 2 + <1 2 , 

eH, = e . qC ° SV + Q(e 2 ). (7) 

^l + ip + Ikf + q 2 

In the main approximation, the equations of motion are 

dip v /l + (p + Ik) 2 + q 2 
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&H k(p + Ik) 



dl ~ y/l + (p + Ik) 2 + q 2 

on q 



(9) 



P 



dq " y/1 + (p + Ik) 2 + q 2 

on p + ik 




dp ' ^/l + ( p + Ik) 2 + q 2 ' 



In this system, variable ip is fast (its time derivative is a value of order 
e), and the other variables are slow (their time derivatives are of order 1). 
Thus, one can average over fast phase ip and obtain the averaged system. 
Motion in this system is just the Larmor rotation in the uniform magnetic 
field B . The averaging, however, does not describe the motion adequately 
near the resonance = 0. At the resonance, projection of particle's velocity 
onto the gi-axis equals the phase velocity of the wave. Resonance condition 
dHo/dl = defines a resonance surface in the (p, q, 7)-space: 



We denote the value of / on this surface as I r . This is a function of variables 
p,q: 



Intersection of surface fflUj) and isoenergetic surface Hq = h defines the res- 
onance curve. Its projection onto the (p, g)-plane is a hyperbola given by 
equation 



where we introduced the dimensionless phase velocity of the wave = u/k. 
Note that is always smaller than 1. 

Variable / is an integral of motion of the averaged system (see (jHJ)) and 
thus an adiabatic invariant of the exact system. Far from the resonance 
the value of I is preserved with the accuracy of order e on time intervals of 
order 1/e (see, e.g., [3]). The adiabatic invariance of / breaks down near the 
resonance, where the averaging does not work properly. In a neighborhood 
of the resonance phenomena of capture and scattering can take place. We 
study dynamics of the system near the resonance in the following sections. 




(10) 




(11) 




(12) 
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3 Motion near the resonance 

To study the system near the resonance, we apply the approach formulated 
in [27] (see also [2H])- Close to the surface ({TO]) the Hamiltonian can be 
expanded into series in (/ — I r ): 

H = A(p, q)+ l -g{v, g)(/-/r(p, g)) 2 +^i|/=/ r +0(|/-/ r | 3 )+0(e(/-/ r (p, g)) 2 ). 

(13) 

Here A(p,q) = H \j = i r is the unperturbed Hamiltonian H Q restricted onto 
the resonant surface (TTUl) . Function g(p,q) in ( ITBl is d 2 H /dI 2 restricted 
onto the resonant surface. It is straightforward to find 

,(l-^) 3/2 



A(p,g)=p^ + J(l + g 2 )(l-^), ^( p ,g) = fc 2 ^-=^. (14) 



Introduce notation d = gy 1 — t^/ a/1 + g 2 . Then 

e-ffi | i=i r — ed cos v?. (15) 

Now we introduce new momentum K = I — J r (p, g) + 0(e). To this end, we 
make a canonical transformation of variables (I,ip,p,q) i— >■ (K,ip,p,q) with 
generating function Wi = pe~ x q + (K + I r (p,q))ip. Omitting bars, we find 
for the Hamiltonian in the new variables (we neglect terms of higher orders) : 

W = A(p, q) + ^g(p, q)K 2 + ed cos tp + eb(p, q)ip, (16) 

where 

6 (p, g = {I r , A} = ——-—— = 17 

<9g dp dp <9g fc^i + g 2 A _ ^2 

Introduce P = K/yfe, 9 = ty/e, and rescaled Hamiltonian F = H/e. The 
rescaled system is Hamiltonian, and pairs of canonically conjugate variables 
are (P,tp) and (p,e~ 3 / 2 q). In the main approximation, the Hamiltonian is 

F = e- 1 A(p,q) + F (P,v,P,q), F = l -gP 2 + dcosy? + hp, (18) 
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Figure 1: Schematic view of the phase portrait of the system (fl8|) for d/b > 1 
(a) and for d/b < 1 (b) 



and the equations of motion are 



P 



-ve 



P' 



dA 
dq ' 
_dFo 

dip 



q 
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dA 

dp ' 

dP 



(19) 



where prime denotes derivative over 9. One can see that variables (P, <p) are 
fast, and variables (p, q) are slow. Thus, as the first step to study this system, 
one can consider variation of the fast variables at fixed values of p and q. 
Dynamics of the fast variables is defined by Hamiltonian Fq, which contains 
q as a parameter. This is a Hamiltonian of a pendulum under the action of 
external torque. Consider phase portrait of this system (fast subsystem). If 
d/b > 1, there is a separatrix on the phase portrait, and, correspondingly, 
domain of oscillatory motion (see Fig.l). In the opposite case d/b < 1, there 
is no separatrix. Note that ratio d/b = k(l — vV) is independent of q. In 

original dimensional units d/b = (B w /Bq)(1 — v^/c 2 ), where = u/k. We 
consider first the case d/b > 1. 

It is straightforward to obtain for the area S inside the separatrix on the 
plane (ip, P) 



S 



<Pm / d 

2- 
9 



cosy? s + -(p s - cos(p -w 

d d 



1/2 



dip, 



(20) 
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Figure 2: Left panel: trajectory of the particle in momentum space (pi,P2), 
Right panel: the same situation as in left panel with a longer time range. 

where ip s = arcsin(6/<i) and ip m is the root of equation cosy? + (b/d)cp = 
cos(p s + (b/d)(p s different from (p s . One can see that both integration limits 
and the expression in square brackets in fj20|) do not depend on q. Thus, 



where A is a constant independent of q (and p). 

Now take into account slow variation of p and q according to the first 
two equations in (1T0]) . It follows from the expression (TT4"]) for A(p,q) that if 

> 0, variable q grows with time. This means that area S also grows in the 
process of motion. Therefore, phase points on the phase portrait of the fast 
subsystem can cross the separatrix and enter the domain of oscillations. This 
corresponds to a capture into resonance. The area J encircled by a trajectory 
in the domain of oscillations on this portrait is an adiabatic invariant (it is 
called the inner adiabatic invariant). Thus, as S monotonously grows with 
time, a captured particle goes deeper and deeper inside the oscillation domain 
and cannot leave it. This means that a particle captured into the resonance 
is captured forever. 




(21) 
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Motion of a captured particle can be described as follows. In the main ap- 
proximation, it moves with the resonant flow defined by Hamiltonian A(p, q). 
The corresponding equations of motion are: 

<9A 
oq 




q = e— = EVcj,. (22) 



It means that a captured particle moves in gx-direction with the wave at 
a speed of the wave's phase velocity. It follows from fl2]) and the fact that 
V2 = that P2 ~ —q, and, hence P2 ~ —ev^t. Therefore pVcomponent of 
the particle's momentum varies (on average) linearly in time: 

p 2 ~ --Bov^t ~ mciojjjv^t/c) 

Thus, the particle is accelerated along the wave front. This acceleration 
is called surfatron one. To find variation of p\ in this motion, we use that 
Pi = p+Ik and expression ( TTTj) for / on the resonant surface. Thus we obtain 

Pi = v+yj(l + q 2 )/(l-vl) and 

n = ^ / . (23) 

where we used the second equation in (1221) . At g ^> 1 we find that p\ grows 
with time as ev^t/ yjl — In dimensional variables, we find that 

Pi ~ l=== t ~ mc / ^={u) L v <t> t/c). 

For the energy of a captured particle = i/l + p 2 + P2 we find that it 
also grows linearly with time at large enough values of q. Namely, we 

have E ~ ev^tj yjl — v 2 and, in dimensional variables, E = mc 2 E ~ 

meruit/ yl - ^/c 2 - 

The captured particle also oscillates in the potential well of the wave. 
These oscillations correspond to motion in the oscillatory domain in Fig. 1. 
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Figure 3: Particle's energy as a function of time. The inner panel shows the 
short time interval before capture. 

One can evaluate the amplitude and the frequency of the oscillations using 
conservation of the inner adiabatic invariant J. For a captured particle J 
equals the area inside the separatrix at the time when the particle crossed the 
separatrix. It follows from the expression for F in (|T5|) that for a captured 
particle 

2 3/2 , ^ r V 2 rz 

J = — yjT^X f"-^-W^)^ (24) 

where fo = Fq/oI does not depend on tp; <pi and (f2 are the minimal and the 
maximal values of (p on a trajectory with fixed values of Fo and q. Thus, 
growth of q results in decreasing of the amplitude of the (^-oscillations. When 
the amplitude of these oscillations is sufficiently small, one can expand the 
Hamiltonian F to obtain 2(F — F ) « gP 2 + d\ cos ip \ ■ (cp — ipo) 2 , where F 
and (po are values of Fo and ip at the bottom of the potential well inside the 
separatrix (see Fig. la). The frequency of oscillations (in terms of rescaled 
time 6) is approximately y~gd\ cos<^o|; an d in this approximation 2irJ ~ 
(F — F )/y/gd. When the particle is captured, J = Jo- Using conservation 
of J along the trajectory of the captured particle, one obtains the scalings 
Atp ~ q- l '\k 2 - w 2 ) 1 / 4 and AP ~ q x ^{k 2 - w 2 )- 1/4 , where Acp and AP 
are amplitudes of yj-oscillations and P-oscillations accordingly. Thus, the 
amplitude of oscillations in q decreases with time proportionally to t -1//4 , 
while amplitude of oscillations in p 1 grows with time proportionally to t l / A . 
Accordingly, amplitude of oscillations in p 2 decreases with time as t~ l / A . In 
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Figure 4: Oscillation frequency of (p as a function of time for two parti- 
cles with different values of initial momentum. Grey dashed lines show the 
theoretical dependence loq ~ t^ 1 / 2 



dimensional variables we have Aqi ~ e(v t f ) i)~ 1 l i and Api ~ y^f/u^t) -1 / 4 . The 
frequency of these oscillations is u>q ~ \fegd (we recall that we made time 
rescaling to obtain ( !T8|) ). Thus, 



uonkil-vftjj^. (25) 

At g > 1 we find that this frequency decreases as >Js/q. In dimensional 
variables we find that co ~ (et)~ l l 2 . 

Here are more complete formulas in dimensional variables for the ampli- 
tudes of the oscillations at large enough values of qi, such that we can put 
qi = v^t. One obtains: 

2\ 1/4 

(c^V/c)- 174 , 



A Ql ~ ec/u Ly /B /B w (l - V -^\ 

Capture into the resonance is a probabilistic phenomenon (see, e.g., [2j 
EH]). Its probability is a small value of order yfe. However, the geometry 
of the system makes the particle pass through the resonance repeatedly, at 
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each Larmor turn. The probability of capture after ~ e -1 / 2 Larmor turns is 
a value of order one (provided that d/b > 1). 

For comparison with theoretical results we present the numerical solution 
of the system with Hamiltonian (j3]), parameters e = 0.1, u/(kc) = 0.25 and 
initial value of the p\ = 0.1. The particle trajectory in momentum space 
(p!,p 2 ) is shown in Fig. [2J At the first stage of modelling the particle 
rotates in the constant background magnetic field (Larmor rotation). This 
motion is slightly perturbed by influence of the wave: the Larmor circles in 
(pi,P2) plane are "scattered". Then after certain time interval the particle is 
captured by the wave and the magnitude of momentum p2 grows with time 
while momentum pi oscillates around the resonant value (it is also increasing, 
yet much slower). To compare the scale of growth of momenta p\ and P2 we 
plot the same picture on a longer time range (see Fig. [2]). The relation 
between p\ and P2 after initial time interval is 1/40 ~ eu/(kc), in agreement 
with the theory (see equation f )23|) ). The particle energy E = a/1 + p\ + p\ is 
shown in Fig. [3j The energy is almost constant before capture (if we neglect 
small scatterings due to wave impacts) and after the capture the energy 
grows linearly with time. In addition we examine the theoretical equation 
for frequency of oscillation of the captured particle - equation ( 125]) . For this 
purpose we plot the oscillation frequency of (p around the null value - Fig. HJ 

4 Scatterings on the resonance 

Capture into the resonance is impossible in the case d/b < 1, when there 
is no oscillatory domain on the phase portrait of the pendulum-like system 
(see Fig. lb). However, in this case the particle energy also changes at the 
resonance crossing. This happens due to scatterings on the resonance. If 
d/b < 1, the average value of the jump in the energy is zero (see [27]), but 
the dispersion is non-zero, and thus diffusive variation of the particle energy 
may be possible. Here we study this topic more attentively. 

When the particle is far from the resonance, its energy is approximately 
constant, because the impact of the wave can be averaged. Thus, to study 
variation of the particle energy we find its time derivative according to equa- 
tions of motion (1191) and integrate it near the resonance. We have 

E = yjl+pl + q 2 = y/l + ip + kiy + q 2 . (26) 
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Using (jSJ) we can write 



E 



d 
df 



(27) 



From ( TIB"]) and ( fl"5]) we find that near the resonance 



E = eoodsimp, 



(28) 



To integrate (1281) we change the integration variable from time t to phase 
<p according to p> = ^dE /dP. Thus we find for variation (jump) of the 
particle energy on one resonance crossing 



where (p* is the wave's phase at the resonance crossing, and q is taken also at 
the crossing of the unperturbed trajectory with energy E and the resonant 
surface. The value of y?* strongly depends on initial conditions and should 
be treated as random. Therefore, change in the particle's energy on the 
resonance is also a random variable. If there is no separatrix on the phase 
portrait in Fig. 1 (i.e., if d/b < 1), the average value of this latter random 
value is zero. An important question is whether these values at successive 
crossings are statistically independent. Expressing q in ( 1291) on the resonance 



via E (we use that at the resonance E — y (1 + q 2 )/(l — we find that 

at E ^> 1 the variation of energy at the resonance scattering is a value of 
order 



On the plane {q,Pi) unperturbed motion of the particle is rotation along 
the circle p 2 + q 2 = E 2 — l. Using Hamiltonian equations of the unperturbed 
motion, one immediately finds that the frequency of this rotation is ujq = e/E. 
The resonant curve on the plane {q,Pi) is a branch of hyperbola pf(l — 
V T) ~ Q 2v 4> = v s w ith p\ > 0. At large enough values of E the trajectory 
of the unperturbed motion crosses the resonant curve at two points. It is 
straightforward to find that the time of motion between these two points is 
a value of order E/e. Consider two successive resonance crossings. Let the 
values of (p at the first and the second crossings be ip^ and p^ accordingly. 
To find p^ — ip* one can integrate equation of motion for ip in ([8]). Thus, 




(29) 





(30) 
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Figure 5: The Poincare sections for several cases with different initial energy 
(e = 0.1). Number of resonance crossings in both cases is 10 6 . The energy 
diffusion is clearly seen. 



one obtains <pl — (p\ ~ E/e. A small variation 5ip\ of the phase ip\ 
results in variation of the energy jump 5je ~ y/eE5(f^ at the first resonance 
crossing. The resulting variation 5ip^ of the phase at the second crossing 
can be found as 5j E d(ip^ — tp*)/dE ~ Sj E /£ ~ 8ip*\/E/y/e. Thus, the 
resulting variation in the phase at the second resonance crossing is much 
larger than 5(p* . Therefore, the values of phase at successive resonance 
crossings are statistically independent. Hence, the jumps in the particle's 
energy at the resonance produce diffusive variation of the energy and its 
unlimited stochastic growth. Note, that in [16] the diffusive growth of energy 
was studied in nonrelativistic case. It was found that, unlike in the relativistic 
case, for a nonrelativistic particle the energy diffusion slows down and finally 
stops at large enough energies. 

One can estimate the rate of the energy diffusion as follows. Consider a 
long trajectory that crosses the resonance N ^> 1 times. Introduce new vari- 
able k = \/E. It follows from (1301) that at every resonant crossing k changes 
by a value of order y 7 ^. If successive jumps in energy are not correlated, it 
follows from f )30|) that typical displacement of k after N resonance crossings 
is Ak ~ yeN. Hence, typical value of energy after N jumps is proportional 
to N: 

E ~ eN. (31) 

Time interval between successive jumps is a value of order of the Larmor 
period. Hence, the time interval corresponding to iV resonance crossings is 
t ~ NE/e. Combining this with f l3T|) . we find that the energy typically grows 



14 




200000 400000 600000 800000 1000000 
N 



Figure 6: The energy of particle ensemble (number of particles is 10 3 ) as 
a function of number of resonance crossings in the system with e = 0.1. 
In subpanel: the number of resonance crossings and averaged energy as a 
function of time. 

with time as 

E ~ eVi. (32) 

We also obtain that the number of jumps (resonance crossings) grows with 
time as iV ~ y/i. 

These results on the energy diffusion of a relativistic particle can be ex- 
amined numerically. For this purpose we construct the Poincare section of 
a particle trajectory in (pi,p2) plane. Points on this plane are plotted with 
time period 2tc/u. One can see that the diffusion in (pi,P2) space become 
stronger as the initial energy of particles grows (Fig. [S]). In Fig. [HI we present 
results of numerics illustrating estimates (13T]) and (|32|) . 
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5 Energy loss due to radiation 



The oscillations of the captured particle across the wave front result in energy 
loss due to radiation. The energy loss has the following rate (see [T]): 



d£_ 2e 2 (d 2 q^ 2 



dt 3c 3 V dt 2 
In the dimensionless variables this expression takes the form: 

dE_ 2r co L /dV ' 



dt 3e 3 c \dt 2 



(33) 



(34) 



where we introduced notation r$ = e 2 /mc 2 . For a captured particle moving 
deep inside the domain of oscillations on phase portrait in Fig. lb, we find 

i 2 2 3/4 

^| = -ek-'A^smcoot = ^^^(l - ^) 9 / 4 shiest. (35) 

Here A(p and uq are amplitude and frequency of small oscillations of the 
captured particle, and K\ is a value of order one; Ki depends on the value 
of q at the instance of capture into the resonance (see Section 3). Thus, we 
find 

dt ~ £K2 c k (l + q 2 ) 2{ V<t>) ■ W 

Here K2 is also a value of order one. In the latter expression, the fraction 
containing q reaches its maximum value f max ~ 0.266 at q = a/3/5. Hence, 
loss of energy due to radiation is maximal at this value of q. 

On the other hand, the captured particle is accelerated along the wave 
front, and thus it gains energy at the rate (see Section 3) 



dE 4 



dt 

Comparing expressions (13"6"I) and (I3T|) we find that if 

r . . „.2n5 



^/^l^l. (37) 



K 2 f max ^uj L k\\ - vif < 1, (38) 

cv<p v 



the radiation cannot stop the particle acceleration. In the opposite case 



K 2 f max —co L k 3 (l - v 2 f > 1, (39) 
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the acceleration can be stopped under the additional condition that the cap- 
ture took place at a value of q smaller than the largest of the two roots of 
equation dE + /dt = d£L/di. 

Note that inequality (139]) can be written in dimensional form as 
r UL (B w /B ) 3 (l — v^/c 2 ) 5 > -00. To be valid, this inequality needs either 
B w very large or very small. Neither so large value of magnetic field, nor 
so small value of phase velocity can be found in physically realistic situations, 
and hence fl39l) cannot be valid in such situations. 

6 Conclusions 

In this paper we considered dynamics of a relativistic charged particle in the 
field of an electromagnetic wave in the presence of a background magnetic 
field. We have described the particle capture into resonance with the wave 
and consequent acceleration using approach of the adiabatic theory of mo- 
tion. During the acceleration the particle's momentum in the direction of 
the wave vector pi and along the wave front p 2 change with time linearly 

(p 1 ~ mc J!i^ ==(uj L v^i/c), p2 ~ —mc(uLV<j>i/c) ). As a result the particle 
V 1 ~ 1 W C 

energy E ~ a/ m 2 c 4 + p\ + p\ grows with time as E ~ mc^w^^/c). The 
estimates of energy loss due to radiation of the accelerated particle show 
that it is not sufficient to stop the acceleration and the particle energy grows 
infinitely. If the condition of capture into the resonance is not satisfied (the 
magnitude of the wave is less than a certain value), particle can nevertheless 
gain energy by scatterings on the resonance. In this case E ~ N, where N 
is a number of resonance crossings. 
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